Caloric curve for finite nuclei in relativistic models 
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Abstract 

In this work we calculate the caloric curve (excitation energy per particle as 
a function of temperature) for finite nuclei within the non-linear Walecka model 
{S) ' for different proton fractions. It is shown that the caloric curve is sensitive to the 
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proton fraction. Freeze-out volume effects in the caloric curve are also studied. 



PACS number(s): 21.10.-k, 21.30.-x, 21.65.+f,25.70.-z 

The production of several intermediate mass fragments in a short time scale during 

<■> ! 

heavy ion collisions is known as nuclear multifragmentation. At several hundred MeV/u, 



the multifragment decay follows the formation of an equilibrated projectile remnant. The 
existence of equilibration is consistent with different experimental observations, such as 
the symmetry of the measured rapidity distributions of fragments with Z > 3 ffl. The 
spectator matter has, therefore, been used to investigate a thermally driven liquid-gas 
phase transition [|], 0, [J. One of the evidences of this transition in infinite systems is 
the fact that the heat-capacity exhibits a peak at a given temperature. However in finite 
systems the situation is more complicated [0, 01 H. The caloric equation of state, which 
is given by the excitation energy per nucleon in terms of the thermodynamic temperature 
is an important quantity to be investigated in the search for a signal of a phase transition. 
Nevertheless, it was recently pointed out that the identification of the existence of a phase 
transition cannot be based only on the behaviour of the caloric curve and a more detailed 
knowledge of the thermodynamic phase diagram is also required 0. In particular it was 
shown that the interpretation of the data is sensitive to the use of a variable free volume 
on the calculation. 
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Recently there has been a big development in the description of nuclei and nuclear 
matter in terms of relativistic many-body theory. In particular, the phenomenological 
models developed using the relativistic mean-field theory describe well the ground-state 
of both stable and unstable nuclei [H, B|. These same models are used to describe the 



properties of neutron stars and super-novae pT3| . Therefore, it is important to test these 
models at finite temperature and different densities. In particular, it would be interesting 
to compare the caloric curve obtained within a relativistic Thomas-Fermi calculation with 
the recent experimental data, and verify whether the proton-neutron asymmetry of the 
hot source gives information on the symmetry energy term of these models. 

Within the framework of relativistic models, the liquid-gas phase transition in nuclear 
matter has been investigated at zero and finite temperatures for symmetric and asym- 



metric semi- infinite systems []T0|, [IT], [12], |13|, [LJJ. With the help of the Thomas Fermi 
approximation, we have investigated droplet formation in the liquid-gas phase transition 



in cold |T5], [L6| and hot |17| asymmetric nuclear matters using the non-linear Walecka 
model (NLWM) [EL TB\. As shown in Refs. [|T5|, [T7J, the optimal nuclear size of a droplet in 



a neutron gas is determined by a delicate balance between nuclear Coulomb and surface 
energies. The surface energy favors nuclei with a large number of nucleons A, while the 
nuclear Coulomb self-energy favors small nuclei. 

In the present work we calculate the caloric curve, given by the temperature dependent 
excitation energy per particle, for the nuclei obtained with the approach mentioned above. 
In particular, we will study the influence of the proton-neutron asymmetry. We chose 
the systems 150 Sm and 166 Sm because they lie in the mass and charge range of interest 
for the experiments we are analysing |2|, 21]. In the first two references the caloric 
curve presented is obtained with prefragments in the mass range 50 ~ 100 to 200. In 
the third reference the data were obtained for a compound nucleus of mass ~ 160. We 
took two isotopes with quite different number of neutrons in order to study the effect of 
proton-neutron asymmetry. In the framework of the Thomas-Fermi theory, shell effects 
are washed out. Hence, we are calculating average properties. We expect that the caloric 
curve of a given system may depend quantitatively on the system mass but the qualitative 
features, namely the dependence of proton-neutron asymmetry and the effect of freeze-out 
volume, will be similar. In multifragmentation calculations an input parameter called the 
freeze-out radius is normally used |4[], so that a phase transition at constant volume is 
simulated. We investigate the consequences on the caloric curve when thermalization in 
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a freeze-out volume is imposed in the present framework. 

We start from the Lagrangian density of the relativistic non-linear Walecka model [ITS] , 
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where 0, V M , and N 1 are, respectively, the scalar-isoscalar, vector-isoscalar, vector- 
isovector meson fields and the electromagnetic field; = dpV v — d v V p , B pv = d^b u — 
d^bfj, — g P (bp x by) and F pu = d p A u — d u A^, with the following parameters: the nucleon 



mass M = 938 MeV, the masses of the mesons m s = 492.25 MeV, m v = 795.36 MeV, 
m p = 763.0 MeV, the electromagnetic coupling constant e = \J~^f and the self-interacting 
coupling constants k and A. The set of constants we use is normally identified as NL1 
||, with C 2 = gfM 2 /ml i= s, v, p, where C 2 = 373.176, C 2 = 245.458, C 2 p = 149.67, 
k/M x 10~ 3 = 2g'l x 2.4578 and A x 1CT 3 = -6g* x 3.4334. This parameterization gives 
a good description of the ground-state properties of all stable nuclei. 

The thermodynamic potential is obtained within the Thomas-Fermi approximation. 
After it is minimized in terms of the meson and electromagnetic fields, the following 
coupled differential equations have to be solved: 



V 2 = m 2 + -k<P 2 + -A0 3 - g s p. 



(2) 



VV = m 2 v V - g v p B , 



(3) 




(4) 



V 2 A 



o = -ep p , 



(5) 



where 
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with 

/i±0,p,t) = — 77- w , i = P,n, (6) 

where Ui — fa — Vio are the effective chemical potentials with fa being the chemical 
potentials for particles of type i and 

V p0 = 9vV + y&o + eA , V n0 = g v V - ^-b ; 

e = a/p 2 + M* 2 , M* = M — g s (f) is the effective nucleon mass and T is the temperature. 
Moreover, p B = p p + p n , P3 = P P ~ pn with 

f d3 P 

Pi = 2 J J2^(fi+ ~ i= P, n , ( 7 ) 

and the energy density, obtained from the thermodynamic potential reads: 

£(r) =2^|^ [e(/ i+ + /,_) + V lQ {f l+ - /«_)] 
4 [(V0) 2 - (Wo) 2 - (vm 2 - (VA ) 2 ] 

The coupled differential equations are solved numerically. For more details on the 
analytical and numerical procedure, please refer to |L5|, [I7| . Three kinds of instabilities 
can occur in this system. The condition for mechanical stability requires that (j^f^J > 
, where P is the pressure and Y p = p v j pn is the proton fraction. The condition for 
diffusive stability implies the inequalities \j£r) > and \j£p-) < . Finally, the 

\ cIY p J pt \ " y p J PT 

thermodynamical stability is expressed by C„ = (^)„y > 0> where C v is the specific 
heat and e* = e(T) —s(T = 0) is the excitation energy per particle, with the total energy 
per particle at any temperature given by e(T) = / ^d 3 r, where A = Z + N. The 
two-phase liquid-gas coexistence is governed by the Gibbs condition. 

We have first solved the equations of motion for an infinite system in order to obtain 
appropriate boundary conditions for the program which integrates the set of coupled non- 
linear differential equations (|2|) to (|5|) in the Thomas-Fermi approximation. Once the 
fields are obtained, all thermodynamic quantities of interest can be easily calculated. The 
binding energy per nucleon is f = e(T) - 938. MeV. 
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In table 1, we show the binding energy per nucleon and the excitation energy per 
particle for the g^Snigg, which has a proton fraction equal to 0.41. In table 2, the same 
quantities are shown for the g2 6 Sniio4, with a proton fraction of 0.37. Notice that, indepen- 
dently of the proton fraction, the excitation energy per particle increases with temperature 
in the range of temperatures shown. For higher temperatures we were not able to obtain 
convergence for a droplet of the size considered. 

In tables 3 and 4, we give the binding and the excitation energies per particle when 
a freeze-out volume of respectively 6V and 9V is used with V the volume at T = 0. 
We have considered a freeze-out radius of 2.2A 1 / 3 fm for 6Vo and 2.5A 1 / 3 fm for 9Vo 
with A = 166. In this case the solutions obtained consist of a droplet immersed in a 
gas of evaporated particles, in such a way that they mimic a source of changing mass. 
As temperature increases more particles evaporate, mainly neutrons, and the fraction of 
protons in the droplet increases. This can be seen in tables 3 and 4 where the number 
of particles which remain inside the droplet as well as the droplet proton fraction (Yd) 
are given. We conclude that the larger the freeze-out volume the faster the excitation 
energy increases with temperature and the larger is the proton fraction in the droplet. 
This picture is consistent with the discussion presented in ||. 

The results for the excitation energies shown in all tables are displayed in figure 1. Also 
shown are the experimental data of refs. |1|, 0, and the Fermi-gas law s* = 1/kT 2 , with 
k = 10.0 (thin dashed line) and 13.0 (thin full line). We have considered that the measured 
temperature Tjj 6 l^ (^HeTD)' obtained from the isotope yield ratios 3 He/ 4 He and 6 Li/ 7 Li 
( 3 He/ 4 He and 2 H/ 3 H), satisfy, in the range of densities considered, T exp /T ~ 0.85 and 
have scaled the experimental data accordingly 0. 

We conclude that the excitation energy for 166 Sm (thick long-dashed curve), proton 
fraction 0.37, increases slightly slowlier with temperature than for 150 Sm (thick full curve), 
proton fraction 0.41, although the difference is not large. These two curves are consistent 



with data of |^T[ and a level density parameter A/k, k = 13.0 in the Fermi gas model 
relation. This agrees with the observed value at around 2 MeV excitation energies [f2~T| . 
Experimental results obtained at higher bombarding energies j|J [J give higher excitation 
energies for the same temperatures. It can be seen from figure 1 of ref. H that the 
higher excitation energies correspond to smaller sources. The larger sources with an 
average -^b ounc i > 60 have excitation energies e* < 5 MeV. In the present approach the 
solutions obtained in a fixed volume, correspond to droplets in a gas of free particles. 



5 



These solutions have higher excitation energies than the ones obtained with no a priori 
fixed volume. In average, this situation corresponds to smaller systems at higher energies. 
This could explain the change of slope that is observed in the calculated data both for 
V = 6 Vq and V = 9 V , in such a way that they come closer to the experimental data 
[|l], 0. The same effect was obtained in ref. p3| , where an exact analytical solution of the 
statistical multifragmentation model was found in the thermodynamic limit. For a fixed 
nucleon density, the caloric curve rises more slowly for lower densities and its leveling 
occurs at lower temperatures. The leveling of the caloric curve is associated with the fast 
change of the configurations from a state dominated by one liquid fragment to a gaseous 
multifragment configuration. We can draw a similar conclusion from tables 3, 4 and figure 
1: the leveling of the thick dash-dotted and short-dashed curves occurs faster when the 
droplet (liquid phase) becomes smaller. It would be interesting to study the effect of the 
symmetry energy on the leveling of the caloric curve in the statistical multifragmentation 
model of ref. [ f23fl . 

In summary, we have studied the excitation energies of arising droplets in a vapor sys- 
tem, up to T = 6.5 MeV. The droplets are described in terms of a non-linear Walecka-type 
model within the Thomas-Fermi approximation. We have used the NL1 parameteriza- 
tion, which is known to describe well the ground-state properties of nuclei. The excitation 
energies of droplets either corresponding to 150 Sm or 166 Sm, for temperatures between 3 
and 6.5 MeV, are consistent with the caloric curve in the Fermi gas approximation with 
a level density parameter A/13. This result agrees with experimental data obtained in 
heavy- ion collisions at intermediate energies We have shown that the caloric curve 
is sensitive to the proton fraction and therefore to the symmetry term of the model used. 
Experimentally the dependence on the proton fraction could be studied by comparing 
data obtained from sources with different proton fractions. For the range of temperatures 
studied, the NL1 parameterization of the non- linear Walecka model has shown to be ad- 
equate to describe nuclear properties and therefore it is a good candidate to generate an 
equation of state for astrophysical purposes. 

When a freeze-out radius is imposed, our procedure yields caloric curves which come 
closer to the experimental results obtained in heavy -ion collisions at higher energies []]. [2] . 
In this case we have smaller droplets with a higher proton fraction immersed in a gas of 
particles, mainly neutrons. This could be interpreted as an oversimplified picture of the 
second regime in the statistical model prediction ||22|| , namely the coexistence phase with 
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a multifragment mixture. This interpretation is supported by the results of ref. . 

Although the thermodynamical equilibrium analysis oversimplifies the problem of high 
energy heavy-ion collisions it is useful for providing a concrete description of warm nuclei 
and for showing qualitative features that should be present in more microscopic calcula- 
tions. 
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Table 1. Output results given by the solution of the coupled differential equations for 
different temperatures for g2°Sm 88 (Y p = 0.41). 



T B/A e*(T) 

(MeV) (MeV/A) (MeV/A) 



0. 


-8.2 


0.0 


3. 


-7.6 


0.7 


4. 


-7.1 


1.2 


5. 


-6.4 


1.8 


6. 


-5.4 


2.8 


6.5 


-4.8 


3.4 



Table 2. Output results given by the solution of the coupled differential equations for 
different temperatures for g2 6 Smio4 (Y p = 0.37). 



T B/A e*(T) 

(MeV) (MeV/A) (MeV/A) 






-7.9 


0.0 


2. 


-7.7 


0.3 


2.5 


-7.5 


0.4 


3. 


-7.3 


0.7 


4. 


-6.8 


1.1 


5. 


-6.1 


1.8 


6. 


-5.2 


2.7 


6.5 


-4.8 


3.1 
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Table 3. Output results given by the solution of the coupled differential equations for 
different temperatures for g2 6 Smio4 in a fixed volume, 6Vq. 



T 


B/A 


e*(T) 


A 


Yd 


(MeV) 


(MeV/A) 


(MeV/A) 






0. 


-7.9 


0.0 


166 


0.37 


3. 


-7.1 


0.8 


163 


0.38 


4. 


-6.3 


1.6 


158 


0.39 


5. 


-5.4 


2.5 


154 


0.40 


6. 


-3.9 


4.0 


148 


0.40 


6.5 


-3.0 


4.9 


143 


0.40 



Table 4. The same as figure 3 for a fixed volume, 9V . 



T 


B/A 


e*(T) 


A 


Y d 


(MeV) 


(MeV/A) 


(MeV/A) 






0. 


-7.9 


0.0 


166 


0.37 


3. 


-7.0 


0.9 


161 


0.38 


4. 


-6.2 


1.7 


156 


0.39 


5. 


-5.1 


2.8 


150 


0.41 


6. 


-3.4 


4.5 


142 


0.41 


6.5 


-2.2 


5.3 


134 


0.41 
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e * (MeV) 

Figure 1: The caloric curves are shown for g2 6 Sm 104 (Y p = 0.37 - thick long-dashed line), 
62°Sm88 (Y p = 0.41 - thick full line); at 6 V (thick dash-dotted line) and 9 V (thick short- 
dashed line) fixed volumes for 166 Sm, and for the Fermi gas law |21| (k = 10.0 - thin 
dashed line and k = 13.0 - thin full line). Experimental results from M (squares), M 
(triangles) and pl| (big stars) are also displayed. 
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